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Abstract 

Background: Harmonic Nanoparticles are a new family of exogenous markers for multiphoton imaging exerting 
optical contrast by second harmonic (SH) generation. In this tutorial, we present the application of Hyper-Rayleigh 
Scattering (HRS) for a quantitative assessment of the nonlinear optical properties of these particles and discuss the 
underlying theory and some crucial experimental aspects. 

Methods: The second harmonic properties of BaTiOs, KNbOa, KiTiOP04 (l<rrP), LiNbOs and ZnO nanocrystals (NCs) 
are investigated by HRS measurements after careful preparation and characterization of colloidal suspensions. 

Results: A detailed analysis of the experimental results is presented with emphasis on the theoretical background 
and on the influence of some experimental parameters including the accurate determination of the nanocrystal 
size and concentration. The SH generation efficiency and averaged nonlinear optical coefficients are then derived 
and compared for six different types of NCs. 

Conclusions: After preparation of colloidal NC suspensions and careful examination of their size, concentration and 
possible aggregation state, HRS appears as a valuable tool to quantitatively assess the SH efficiency of 
noncentrosymmetric NCs. All the investigated nanomaterials show high SH conversion efficiencies, demonstrating a 
good potential for bio-labelling applications. 



Background 

Imaging techniques are increasingly accompanied by 
nanotechnological developments aiming at the reliable, 
scalable, and efficient production of nanoparticle-based 
bio-markers with specific physical, optical, and chemical 
properties [1]. The use of multiphoton microscopy com- 
bined with fluorescent labels has undoubtedly led to 
improvements in terms of spatial resolution, depth pene- 
tration in biological tissues, and reduction of photo- 
damage [2,3]. However, these optical labels generally suffer 
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from bleaching and blinking limiting optical contrast and 
observation time. A new generation of nanoparticles based 
on nonlinear optical response is being developed in this 
context to overcome these limitations and complement 
fluorescent labels. Second Harmonic emission from non- 
centrosymmetric nanocrystals (NCs) indeed presents 
appealing features for bio-imaging such as a long-term 
photostability and virtually unlimited wavelength flexibility 
[4-9]. 

In linear optics, we define the macroscopic polarization 
P induced on a medium by an electromagnetic field E 
(characterized by an oscillation frequency co or a wave- 
length X = In c/co) as P oc /£. According to this relation- 
ship, the polarization is proportional to the incident field 
and the linear response of the medium is characterized by 
its susceptibility /. The radiated intensity is proportional 
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to the incident intensity and the frequency remains 
the same. In nonlinear optics, the intensity of the 
incident electric field being much higher, the induced 
polarization is no longer proportional to the incident field 
and additional nonlinear terms must be included, 
cx: + x'^'fi^ + X^^'fi^ + • • • ' where higher-order suscept- 
ibilities / (2), / are tensors referring to the medium char- 
acteristics which determine its nonlinear optical response. 
SH generation is a second-order optical process depending 
on where a fraction of the incident light wave is con- 
verted into radiations with frequency twice the incident 
one. It is worth noting that bulk SH generation is only 
observed in noncentrosymmetric media, limiting the 
choice of harmonic nanoparticles to a specific class of 
materials. 

Optical contrast based on frequency doubling process 
can bring several benefits for bio-imaging. First of all, 
involving exclusively virtual electronic states, it does not 
imply energy absorption preventing bleaching usually 
observed with fluorescent and luminescent nanoparticles. 
Therefore, it is possible to achieve observations over long 
time durations with no decrease in signal quality [10]. In 
addition, SH generation is a non-resonant process which 
can occur for any excitation wavelength [10-12]. 
Consequently, excitation wavelengths can be selected in 
spectral ranges where the biological tissue absorption and 
scattering are low so that photo-damage is limited and 
penetration depth increased [13]. A wise choice of the exci- 
tation wavelength also allows avoiding hindrance by sample 
autofluorescence increasing image contrast. In terms of 
emitted intensity, the SH signal scattered by NCs is related 
to the square intensity of incident power. Detection and 
imaging of such individual bio-markers generally requires 
high peak-power excitation beams, a condition nowadays 
easy achievable with femtosecond lasers even for very mod- 
erate pulse energies [13,14]. Practically, the preparation of 
SH nanocrystals as bio-labels involves a good knowledge in 
chemical synthesis and functionalization. Nanocrystals in 
the 10-100 nm size range with good crystallinity [15] and 
low cytotoxicity [10,16,17] are typically required as well as 
the development of conjugation techniques [18-20] for an 
efficient interaction between the NCs and the biomolecules 
of interest. Finally, the selection of optical labels with high 
SH efficiency requires the development of experimental 
techniques to quantitatively assess the nonlinear optical 
response of nanocrystals. SH generation is indeed well 
known in bulk crystals and can also be easily characterized 
in thin films and agglomerated powders. However, charac- 
terization techniques for nonlinear nanocrystal ensembles 
are sorely lacking. In this paper, we propose to apply the 
Hyper-Rayleigh Scattering (HRS) technique, originally 
developed for the investigation of harmonic light scattering 
in molecular solutions, in order to characterize the SH 
efficiency of solvent-dispersed non-centrosymmetric NCs. 



The theory and first experimental detection of 
Hyper-Rayleigh Scattering were first exposed in 1965 
[21,22]. In a series of experiments, HRS was later used 
to determine the second-order nonlinear optical proper- 
ties of molecules including species with no permanent 
dipole moments [23,24]. The HRS intensity is defined as 
the incoherent sum of SH signals scattered by different 
sources excited by an incident field of intensity 1^: 
^HRS = GY.ii^iP> < Pi >Vi > where the subscript i 
refers to the different species in solution. < j8 > stands 
for y< /J2 ^ which is the average amplitude of the 
second-order molecular hyperpolarizability and brackets 
indicate isotropic orientational averaging. This quantity 
describes the molecular nonlinear optical response con- 
sidering that local microscopic fields Eio^. induce dipole 
moments pim that can be written as p2w = fiE^c where P 
is here a third-rank tensor. AT, is the species concentra- 
tion, F a local field factor depending on the surrounding 
medium refractive index, and G an experimental 
proportionality constant. From the above HRS intensity 
equation, the square dependence on the incident inten- 
sity and the linear behaviour with the concentration 
species can be noticed. Besides, quantitative assessment 
of the hyperpolarizability < jS > of nonlinear species 
with the external reference method is specifically based 
on this linear trend [23,24]. By varying the molecule 
concentration (NO and measuring the corresponding 
HRS intensity (Ihrs)- one can obtain the slope a 
defined as a = GP < fi-^ > from the /hrs(M) curve 
a = GF < P'^ > I^. Hyperpolarizability is then derived by 
comparison with the experimental slope measured from 
a reference solution containing molecules of known 
hyperpolarizability. The above procedure assumes that 
the parameter G remains unchanged for defined experi- 
mental conditions whereas local field factors can be cal- 
culated from solvent indices. 

Methods 

Second harmonic NC characterization from HRS 

HRS has been recently extended to larger objects such as 
molecular assemblies or nanometric- to microscopic-scale 
particles [25-31]. A delicate issue about nanocrystal char- 
acterization is that the microscopic and macroscopic 
description of the SH optical process seems inappropriate. 
As previously mentioned, HRS is usually employed to 
determine the molecular hyperpolarizability < >. 
On the other hand, SH properties of bulk crystals are 
commonly associated with the macroscopic second order 
nonlinear susceptibility / (note that / 1^) = 2d, where d 
is the nonlinear optical coefficient in the convention I 
used here [32]). An effective hyperpolarizability can be 
therefore conveniently introduced with the purpose of 
relating microscopic to macroscopic parameters allowing 
a straightforward comparison with the bulk nonlinear 
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optical properties. Assuming that NCs are large enough to 
neglect surface contribution of the scattered SH signal, the 
effective hyperpolarizability < fine > has the same proper- 
ties as X and is related to the averaged SH coefficient 
<d>hY the mean particle volume Vnc [30,33] : 



>=< d > V„ 



(1) 



It should be mentioned here that, contrary to molecular 
hyperpolarizability, < pnc > is not defined according to 
local microscopic fields. An internal field factor T must 
however be introduced to better estimate the incident 
optical macroscopic field in the nanocrystal. In the case of 
colloidal suspensions of nanocrystals, there is indeed a 
refractive index difference between the solvent and the 
investigated nanomaterial. According to this definition, 
the HRS intensity of a given suspension containing NCs 
(subscript "nc") dispersed in a solvent (subscript "s") can 
be expressed as: 



lHRS=G{NsFs<fi^>+N„cT„c<:Pl>)ll 



(2) 



The external reference method can again be applied to 
determine < fine > provided that the suspension concen- 
tration is known. In addition, quantitative assessment of 
the < d > parameter also assumes that the nanocrystal 
size is independentiy estimated. 

Our HRS experimental set-up [34] is based on a 
Q-switched Nd:YAG laser (Wedge HB, Bright solutions, 
1 mJ/1 ns). The vertically polarized beam is focused by 
a 20 cm focal length lens into a glass cuvette containing 
the NCs suspension. Scattered SH signal is collected at 
90° from the laser beam axis with a 5 cm focal length 
lens and measured by a photomultipher tube. A colored 
glass filter (short pass) and an interferometric filter (532 
nm) are placed before the photomultipher to remove 
unwanted signals. The unpolarized HRS signal is finally 
processed by a boxcar amplifier after integration and 
averaging over typically 1000 laser shots. 

Preparation of NC suspensions 

Different types of particles are investigated in this work: 
agglomerated BaTiOs and KNbOs nanopowders obtained 
by ball grinding of bulk crystals (kindly supplied by FEE 
GmbH, Idar-Oberstein, Germany), KTP nanocrystals 
extracted from raw powder at the end of the flux-growth 
process (kindly provided by Cristal Laser, Messein, France) 
[35]. LiNbOs nanocrystals synthesised by a hydrolysis 
method [36] supplied by the SRSMC laboratory (Nancy, 
France), and, two commercial batches of ZnO with 
different average nominal sizes of 20 nm and 90-200 nm 
purchased from NanoAmor Inc. (Houston, USA). 

Colloidal suspensions with low size dispersion and good 
stability are obtained following a specified protocol 
[10,37]. Dry nanopowders are first dispersed in water or 



ethanol at 0.5 mg/mL concentration and sonicated 
(Vibra Cell 75043, Bioblock Scientific) during 25 min in 
pulsed mode (1 s on, 4 s off). Large particles and 
remaining aggregates are then excluded by leaving the 
solution settling down during a period of 1 to 7 days 
depending on the NCs type. As-obtained supernatants 
have a typical concentration of 0.1 mg/mL. Size and size 
dispersion are then characterized by Dynamic Light 
Scattering (Malvern ZetasizerNanoZS). Concentration 
of solvent-dispersed NCs is measured after evaporation 
of a precise volume of the obtained suspension and 
weighing the residual. 

Results 

Hyper-Rayleigh Scattering Data Analysis 

AHRS intensity measurements are performed according to 
the concentration of NCs and, on the other hand, of para- 
nitroaniline (pNA) diluted in methanol which is used here 
as the reference molecular solution. After linear fitting of 
the experimental coefficients oimc and oipNA (that relate 
Ihrs to the nanocrystal and pNA concentrations, 
respectively), the effective hyperpolarizability < fine > is 
calculated as [30,33]: 



OCncFpNA 
OipNATne 



pNA' 



(3) 



FpNA, the Lorentz-Lorenz local field factor correction, is 

Keth + 2 



given by R 



'pNA 



with nmeth the refractive 



index of methanol. The internal field factor is expressed as 

/ 3n^ \® 

= I * — I where «s and ^nc are the refractive 

indices (considered here without index dispersion i.e., 
n{co) ~ n{2(o)) of solvent and NCs, respectively. The 
reference value of the pNA molecular hyperpolarizability 
is /?33 = 25.9 X 10^^''esM[38] from which the orientation- 



ally averaged hyperpolarisability corresponds to y 

in our experimental configuration (vertical input polariza- 
tion and no analyser). 

To allow a comparison with bulk crystal literature 
values, the averaged nonlinear coefficient <d> is then 
obtained according to equation (1) by inserting the 
volume V„c. This procedure applied to the 6 available 
types of NCs, leads to the experimental determination 
of the effective hyperpolarizabilities and averaged 
nonlinear coefficients as shown in table 1. The validity 
of the presented protocol is discussed hereafter with a 
special focus on the possible optical artefacts associated 
with the SH measurements. Experimental precautions 
and recommendations are also introduced with regard 
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Table 1 Results 





Mean DLS size [nm] 


< p„c >[10'^'' esu] 


<d„ono> [pm/V] 


<c'poiy> [pm/V] 


<clbuik>* [pm/V] 


BaTiOs 


91 


5.16 


5.5 




14.1 - [45] 


KNbOj 


128 


21.8 


8.3 


5.6 


1 5.3 - [46] 


KTP 


237 


33.5 


2.8 


1.4 


7.6 - [46] 


LiNbOj 


125 


16.3 


74 


4.9 


1 7.3 - [47] 


ZnO 90-200 


132 


7.12 


2.5 


1.2 


2.8 - [48] 


ZnO 20 


50 


0.33 


2.1 




2.8 - [48] 



Mean NC size considering DLS distribution by number, effective liyperpolarizability < ySnc >> averaged nonlinear optical coefficient derived from monodisperse 
(< dmono >' ^""^ polydisperse « dpoly » suspensions and comparison with literature values (< dimlk >' different SH nanocrystals. Error bars for 

^ f^nc ^ < dmono ^ estimated at ± 15% in terms of experimental reproducibility. 
*<dS = 6/35 d33^ + 32/105 dji.djj + 92/105 dj,^ for BaTiOj (4 mm) and ZnO (6 mm) 

<dS = 6/35 d33^ + 38/105 (djj^ + dj,^) + 16/105 (djj.djj + dji.djj + djj.dsj) for KTP (mm2) and KNbO, {mm2) 
<d^> = 6/35 d33^ + 92/105 dsi^ + 8/21 d22^ + 32/105 dsi.djs for LiNbOs (3 m) 



to the influence of the size, size distribution and aggre- 
gation state of the solvent-dispersed NCs. 

Optical considerations for a correct HRS measurement 

Among the optical considerations and possible artefacts, 
let us first point out that the term HRS is specifically used 
here because scattering at 2(0 is detected with the classical 
HRS experimental configuration. According to the recent 
review of Roke et al. [31] Second Harmonic Scattering 
would be more appropriate for -100 nm non-centrosym- 
metric NCs. If an incoherent process is underlying in eq.2 
(Ihrs increases linearly with the NC concentration), the 
SH intensity also depends quadratically on the NC volume 
as it can be seen after combination of eq.l with eq.2. The 
number of induced nonlinear dipoles (at the scale of the 
crystalline unit cell) being proportional to V^c, a coherent 
process thus takes place within each NC between the var- 
ious radiated SH fields. This bulk response is equivalent to 
the one observed in single crystals except that phase 
matching constraints are absent because the NC size is 
well below the coherent length [5]. Coherent effects are 
thus to be considered within each NC whereas the inco- 
herent contribution originates from the linear concentra- 
tion dependence. The Rayleigh regime assumption, at the 
very basis of HRS theory could be verified for most NC 
suspensions. The scatterers are assumed to be point-like, i. 
e. much smaller than the excitation wavelength so that 
light is scattered isotropically. An experimental verification 
can be carried out by detecting the normalized HRS 
intensity of a NC suspension at different collection angles 
spanning from 90° to 180° with respect to the incident 
beam direction. Figure 1 shows the curves obtained for 
ZnO and KTP. The ZnO suspension (130 nm average 
particle diameter) displays an almost flat angular response 
leading to conclude on the validity of the Rayleigh regime. 
The KTP suspension, on the other hand, containing larger 
NCs of around 200 nm presents higher forward scattering 
of Ihrs evidencing the limit of the Rayleigh regime 
assumption in this size range. 



In regard to the possible optical artefacts, only the non- 
linear elastic scattering is to be measured and a mono- 
chromator is necessary to spectrally resolve the HRS 
signal and to detect the presence of any possible inelastic 
process. This basic verification has been done for all the 
NCs under consideration. We only observed second 
harmonic scattering and no fluorescence [39]. Light 
absorption and scattering can also lead to alteration of the 
HRS signal. For instance, light absorption at the working 
wavelengths (1064 and 532 nm) tends to reduce the SH 
intensity when the species concentration is too high. This 
results in a deviation from the linear HRS response and a 
Beer-Lambert correction factor is to be added to correcdy 
estimate the expected proportionality between the HRS 
intensity and the concentration [40]. Similarly, a deviation 
has also been observed with non-absorbing NCs (see 
Figure 2) so that an exponential correction was introduced 
to account for the very high multiscattering processes of 
the incident and SH radiations. Note that further dilution 
of the initial supernatants readily allows the recovery of a 
linear dependence on the NC concentration. 

Influence of the NC size, concentration and possible 
aggregation state 

Size characterization of NCs is of paramount importance 
for the correct determination of SH efficiencies. As 
previously described, the averaged nonlinear optical coeffi- 
cient <d> is retrieved from < p„c > and the mean NC 
volume. Combination of eq.l with eq.2 shows a depen- 
dence of HRS intensity with the squared NCs volume V^^, 
emphasizing the requirement of a precise size characteri- 
zation. In addition, assessment of the suspension concen- 
tration also relies on the NC size. As detailed hereafter, a 
mass concentration is indeed first measured and the unit 
particle concentration then estimated according to the NC 
volume. 

Size characterization can be performed using Dynamic 
Light scattering (DLS), a technique based on the intensity 
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90 100 110 120 130 140 150 160 170 180 

detection angle 0 



Figure 1 Rayleigh. Variation of the HRS intensity of ZnO and KTP suspensions for a collection angle spanning from 90° to 180° with respect to 
the incident beam direction. Strong fluctuations observed between 140° and 160° are attributed to a non-significant interaction of the incident 
laser beam with optical components the rotating detection arm. 
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fluctuations in the scattered light. These fluctuations are 
caused by Brownian motion and can be related to the par- 
ticle size. This method provides a particle size distribution 
in intensity since the detected signal is proportional to a 
linear elastic scattering. A main issue may arise from size 
polydispersity because large particles can readily hinder 
the lower intensity response of smaller scatterers. In 
Figure 3, two different LiNbOs suspensions are character- 
ized by DLS. Sample A was prepared as described in 
"methods" whereas the sedimentation step was not 
completed for sample B. A comparison is made with the 
distributions obtained from a particle-by-particle investi- 
gation technique called Nanoparticle Tracking Analysis 
(NT A) [41]. The good agreement between both results for 
sample A demonstrates the reliability of DLS measure- 
ments in that case. On the contrary, the larger dispersion 
associated with suspension B results in a less consistent 
correlation. This again points out the importance of 
paying a particular attention to suspension preparation. 



to 

Q 




100 

Size [nm] 



Q 




100 

Size [nm] 

Figure 3 Size. DLS size distributions by intensity and number for 
two different suspensions A (a) and B (b) of LiNb03 NCs. 
Comparison is made with the NTA distribution by number. 



Other critical issues complicating the interpretation 
of DLS results are related to the nanoparticles shape 
and the presence of aggregates. DLS indeed detects the 
hydrodynamic volume of objects and uses spherical 
models making the size characterization less relevant 
for non-spherical particles. Complementary observa- 
tions of NCs by electron microscopy may help to 
obtain more information about the NC shape polydis- 
persity. In particular, the discrimination between even- 
tual aggregates and large individual NCs is not possible 
with DLS. As an exemplary result, ageing of the initial 
water-based solutions may result in a slow aggregation 
if electro-steric repulsions are too weak after the initial 
preparation of NC suspensions. From our experience, 
the long term stability of water-dispersed NCs is parti- 
cularly difficult to achieve with BaTiOs, ZnO and KTP 
NCs and the use of ethanol as a solvent is generally 
preferred. 

Influence of the aggregation state on the HRS signal is 
specifically investigated in Figure 4. Simultaneous DLS 
and HRS measurements for ZnO NCs suspended in 
ethanol have been carried out for several days. The size 
increase in DLS measurements is typical of a sponta- 
neous aggregation process (Figure 4a) whereas the 
almost constant HRS signal attests of the low effect of 
aggregation on the SH intensity (Figure 4b). This experi- 
ment unambiguously demonstrates that the HRS signal 
only depends on the individual nanocrystal size. Even 
after aggregation, the incoherent contribution of each 
individual nanocrystal is preserved indicating the 
absence of specific orientations within aggregates. 
On the contrary, the individual orientation and dynamics 
of spontaneous aggregation were observed to be well dif- 
ferent for BaTiOs NCs previously dispersed in a nonpolar 
fluid, lilce heptane [42,43]. 

Finally, accuracy of concentration measurements 
should also be clarified. In our experience, the weighting 
method used to determine the mass concentration is 
satisfactory provided that the evaporated volume is suffi- 
cient. Typically, a 30 mL suspension is evaporated in 
several successive steps in light containers of 7 g. The 
resulting relative uncertainties in the mass concentration 
are between 1 and 5%. N„c, the unit particle concentra- 
tion, is then calculated from the mass concentration Cm 

Cm 

with Nnc = where P is the material density and 

pVnc 

V„c the mean NC volume. As individual NCs are moni- 
tored with the NTA analysis (for which the preparation 
of very diluted suspensions is a prerequisite), compari- 
son between the two approaches is straightforward. 
Deviations by a factor of 1.5 were noticed in the best 
cases, which is acceptable considering the inherent size 
polydispersity that affects the accuracy of the weighing 
procedure. 
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Figure 4 Aggregation. Effect of the aggregation state of ZnO NCs on the measured HRS intensity: a) time-dependent DLS measurements 
showing the apparent aggregate size, b) HRS intensity evolution with aggregates size. 



Influence of the NC size distribution 

Another point of interest is that the sizes used for deri- 
vation of the nonUnear optical parameters are actually 
the mean values of distributions. In other words, HRS 
intensity is considered as a nonlinear elastic scattering 
from exactly monodisperse samples. Because of the 
already mentioned size polydispersity, the HRS signal 
(that depends quadratically on V„c) indicates that a 
mean value is not commensurate with the real scattered 
SH intensity. In order to evaluate the influence of the 
actual size distribution width, we have compared the 
HRS intensities of a virtual (or ideal) monodisperse sus- 
pension to an actual polydisperse one. To begin, we 
only take the volume into consideration. We also 
assume that both distributions of spherical NCs (with 
diameter D) have the same mean size, Dm, but a non- 
zero standard deviation for the actual suspension. To 
that end, DLS distributions by number are fitted by 
using a normal probability density function (p (D)). 
Ratios of the HRS intensities between polydisperse and 
monodisperse samples can be then calculated as: 

Jpoly 1 



In 



D6 



/ D^p{D)dD. In the case of the LiNbOg 



suspension A of Figure 3a, the calculated HRS intensity 
(with a 125 nm mean diameter and a 25 nm standard 
deviation) is for instance twice higher than a monodis- 
perse suspension containing exclusively 125 nm NCs. 

The statistical approach being introduced, our final 
aim is to evaluate the effect of the NC size distribution 
on the derived averaged nonlinear optical coefficient 
< d >■ Concentration (Nnc) and volume (Vnc) are thus 
to be considered. Combination of eq.l with eq.2 allows 
to write < >a ha>/[Nnc ■ V^). As discussed above for 
the concentration determination, Nnc is inversely pro- 
portional to the NC volume so that the effect of a 



size-dispersion on the calculated concentration can be 
obtained from — = — - — . Taking into 

Nmono f D3p{D)dD ^ 

account the squared volume, the final effect of a size 
dispersion for spherical NCs with diameter D on the 
averaged < d > coefficient can be expressed as: 

< dpoly > p (D) dP 

<dmono> V "'l°'p{D)dD' 

Interestingly, the whole effect of the NC-size dispersion 
for the LiNbOs A suspension is relatively weak. A factor 
of about 0.8 is indeed found between the "polydisperse" 
and the "monodisperse" averaged coefficient showing that 

< > is generally overestimated when the size distribution 
is ignored. 

Finally and after consideration of the previous recom- 
mendations and precautions, our experimental results for 
the 6 different NC type are summarized in table 1. The 
mean NC size (D^) is obtained from a DLS distribution by 
number whereas the effective hyperpolarizabilities < fine > 
are derived from eq.3. The averaged nonlinear optical 
coefficient calculated with monodisperse < dmom > and 
polydisperse < dpoiy > (when available) suspensions are 
then compared with the literature values < di^ik > of the 
different material. Note that literal expressions of 

< dituik > (corresponding to our experimental polarization 
configuration) can be found in footnote of table 1. They 
have been calculated assuming Kleinman conditions and 
the symmetry class of each material [24,44]. 

Discussion 

Among the different experimental optical parameters, let 
us remind that < fi„c > is characteristic of the NC SH 
efficiency. The effective hyperpolarizabilities < fine > 
indeed linearly depend on the NC volume whereas the 
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averaged nonlinear optical coefficient < > is an intensive 
physical property of each material. This distinction is read- 
ily observed with the two ZnO samples that display similar 

< d > coefficients confirming that the bulk contribution is 
predominant in the measured HRS signal. Consistently, 
large-size NCs (ZnO 90-200) have a higher hyperpolariz- 
ability than small NCs (ZnO 20). For KTP NCs, the high- 
est measured effective hyperpolarizability of 33.5 x 10~^* 
esu is to be related to the relative large NC volume and 
not to especially efficient intrinsic properties. 

In addition, it can be noticed that the < d > coefficients 
for the different materials are otherwise very similar (a few 
pm/V) and that the comparison with literature < dbuik > 
values obtained from bulk crystals is also very consistent. 
Such an agreement is in favour of the HRS technique as a 
quantitative characterization method and also confirms 
the absence of a significant surface contribution in the 
scattered SH intensity. On the other hand, if the apparent 
hierarchy between bulk properties of the different materi- 
als seems preserved, a systematic discrepancy can be 
noticed. Literature bulk crystal and experimental NCs 

< d > coefficients indeed differ by a factor 2 to 3 and sev- 
eral explanations can be suggested. Our quantitative 
assessment is first based on the external reference method 
with pNA molecules for which various different nominal 
hyperpolarizabilities have been published. Similarly, the lit- 
erature nonlinear optical coefficients differ for bulk mate- 
rials. For the NCs, the possible presence of impurities and 
amorphous phases has not been evidenced through X-Ray 
diffraction (data not shown) but slight deviations from the 
chemical stoichiometry could also alter the SH properties. 
Finally, it is our opinion that this systematic shift most 
probably results from the inherent size- and shape- 
distributions together with the likely presence of aggre- 
gates that inevitably leads to an underestimation of the 
nonlinear optical coefficient < d >■ The likelihood of dif- 
ferent ferroelectric domains in the investigated materials 
(except ZnO for which the experimental and literature are 
well consistent) also points out to the necessity to improve 
current strategies for producing NCs with accurate size, 
shape and composition. 

Conclusions 

Application of the HRS technique on NC suspensions 
allows to achieve a reliable assessment of the SH efficien- 
cies of noncentrosymmetric nanomaterials in the context 
of their promising application for bio-imaging as nonlinear 
optical markers. Limitations and experimental caveats of 
the method have been thoroughly discussed, with special 
attention to the crucial size and concentration issues. Our 
conclusions emphasize the importance of the preparation 
of aggregate-free NC suspensions with ideally monodis- 
perse samples in terms of size and shape. These recom- 
mendations should rapidly improve accuracy of these 



quantitative measurements. To date, we can however 
assert that the different investigated NCs all present very 
similar nonlinear optical coefficients < d > and thus com- 
parable SH efficiencies at a given NC size. A pertinent 
selection of a NC type for bio-imaging also weakly 
depends on the nonlinear optical properties. In the future, 
the reliable and scalable chemical synthesis of size- and 
shape-controlled NCs as well as their cytotoxicity and 
biocompatibility are the relevant parameters to consider. 
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